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Abstract 

A combination of the variable-constant and complex coordinate rotation 
methods is used to solve the two-body Schrodinger equation. The latter is 
replaced by a system of linear first-order differential equations, which enables 
one to perform direct calculation of the Jost function for all complex momenta 
of physical interest, including the spectral points corresponding to bound and 
resonance states. Explicit forms of the equations, appropriate for central, 
noncentral, and Coulomb-tailed potentials are given. Within the proposed 
method, the scattering, bound, virtual, and resonance state problems can be 
treated in a unified way. The effectiveness of the method is demonstrated by 
a numerical example. 

PACS numbers: 03.65.Nk, 21.45. +v, 02.70.+d 

I. INTRODUCTION 

The nonrelativistic quantum mechanical two-body problem can be described conveniently 
and nicely in terms of the Jost functions and Jost solutions of the Schrodinger equation. 
When defined for all complex values of the momentum, these functions contain complete 
information about the underlying physical system. In contrast to the conventional treatment 
of the two-body problem, based on the scattering amplitude and physical wave-function, in 
the Jost function approach the bound, virtual, scattering, and resonance states are treated 
on an equal footing and simultaneously. Also the Jost function is a more fundamental 
quantity than the S function since it is free of ambiguities caused by redundant zeros, while 
the S function may have redundant poles 0. 

Despite the fact that the Jost functions are more appropriate, in principle, for studies 
of spectral properties of two-body Hamiltonians, they have not been used in the past in 
practical calculations. Instead they were considered as purely mathematical entities of for- 
mal scattering theory, which can be expressed either via an integral containing the regular 
solution or via a Wronskian of the regular and Jost solutions [0. Thus, to find the Jost 
function based upon such information, the full solution of the problem itself is required. 

In Ref. @, linear first-order differential equations for functions closely related to the Jost 
solutions were proposed. These equations are based on the variable-phase approach || and 
their solution at any fixed value of the radial variable r, provides the Jost function and its 
complex conjugate counterpart, which correspond to the potential truncated at the point r. 
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The proposed method, however, is suitable only for bound and scattering states calculations, 
that is for calculations in the upper half of the complex momentum plane. 

Herein we extend the method of Ref. H onto the unphysical sheet so that the resonance 
state region is included. We do so by a combination of the variable-constant method || 
and the 'complex coordinate rotation' method ||. The latter method, often used in studies 
of atomic and molecular physics, is also known as the 'complex-scaling' or the 'complex- 
dilation' method. It is based on the fact that the complex coordinate transformation, 
f — > fexp(i8) with a sufficiently large 6 > 0, makes the resonance state wave function 
quadratically integrable leaving the energies and widths of the bound and resonance states 
unaffected since the Jost function, and the positions of its zeros in particular, do not depend 
on that vector r. This is true if the potential is an analytic function of the coordinate r and 
vanishes rapidly enough when that coordinate goes to infinity along the turned ray 0. In 
atomic and molecular physics, this property has been used to locate resonances by the same 
variational methods employed for the bound states [0. 

In this paper we make a complex rotation of the radial variable in the differential equa- 
tions proposed in Ref. ||. The new equations that result, facilitate evaluation of the Jost 
function in the entire complex momentum plane and therefore one can locate the resonances 
as the zeros of the Jost function. 

The paper is organized as follows. In Sec. II we derive the rotated equations in the 
simplest case of a central short-range potential. In Sec. Ill the equations are generalised for 
the case of a noncentral potential, while Sec. IV is devoted to consideration of potentials 
having a Coulomb tail. In Sec. V the method is tested by comparing results with the 
properties of an exactly solvable model. Finally, possible further extensions of the method, 
its advantages and drawbacks are discussed in Sec. VI. 



II. CENTRAL SHORT RANGE POTENTIAL 



Consider the problem of two quantal particles interacting via a central potential V(r) 
with the asymptotic properties 

lim r 2 V{r) = , (1) 

r— >0 

and 

lim r 1+£ V(r) = , e > . (2) 

r— >oo 

Any physical solution of the radial Schrodinger equation (h = 2m = 1) 

D t (k,r) $ t (k,r) = V(r) $ t (k,r) , (3) 

where 

D t (k, r) = d 2 r +k 2 - £(£ + l)/r 2 , (4) 

is proportional to the so-called regular solution which is defined by the boundary condition 
at the origin, 
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]3m[$t(k,r)/j t (kr)] = 1 . 



(5) 



Here ji(z) is the Riccati-Bessel function [||. At large distances the regular solution is a 
linear combination of Riccati-Hankel functions hj^\kr) 0, i.e. 



r— >oo 2 



/4 +) (*r) //(*;*)+ /ii _) (Ar) 



(6) 



wherein the r— independent quantity fi(k) is the Jost function whose zeros in the complex 
fc-plane correspond to the bound, (Bek — 0, Imk > ), virtual, (ife/c = 0, im/c < ), and 
resonance, (ifeA; > 0, Imk < ), states of the system. 

Suppose V(r) is an analytic potential @, i.e. it can be defined as an analytic function 
of the complex variable 

r = xexp(z0), x > 0, 9 6 [0, 8 max ], 9 max <ir/2, (7) 

and fulfills the conditions of Eqs. flU) and (@). We consider the complex dilation of the 
regular solution using this rotated coordinate. The corresponding function Qi(k,xe ld ) is 
square integrable over the variable x, even for a resonance state, provided the angle 9 is 
greater than ip which defines the position of the resonance on the lower half of the complex 
fc-plane by 



ko — |&o| exp(— iip). 



(8) 



This is obvious because the product 'fc r' in the asymptotic behaviour ~ exp(i/c r) of the 
solution has a positive imaginary part. 

To obtain the equation for the dilated regular solution, we substitute Eq. (|7]) into Eqs. @ 
and (0) giving 



D e (ke l0 ,x) $ e (k,xe lf) ) = e m V{xe i0 ) $ e (k,xe i(> ) 



(9) 



and 



lim 

x^0 



<5> e (k,xe ie ) /ji{kxe 



1 



(10) 



Furthermore, and since the potential is assumed to be analytic, substitution of Eq. (0) into 
Eq. (P) gives the asymptotic form 



§i(k,xe l 



1 

x— >oo 2 



hf\kxe w ) //(A;*) + h ( f>(kxe w ) f e (k) 



(-)< 



(11) 



Note that the Jost functions in Eqs. (|j) and (|Tl]) are the same because the Jost function 
does not depend on the rotation angle 9. Hence, the positions of its zeros are ^-independent 
as well. 

To solve the differential equation, Eq. (||), we apply the variable constant method ||, 
i.e. we look for a solution in the form 



1 



77l(+)/ 



.(-)/ 



(12) 
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which is an r-dependent combination of two linearly independent solutions of the free equa- 
tion 



D(.{ke l \x) hf\kxe l °) =0 



J9\ 



(13) 



Instead of one unknown function we have introduced two unknown functions which 
necessitates an additional constraint condition. The most convenient is the standard La- 
grange condition || 



hf\kxe l °) d x F^ , {k,x,6) + h [ f > {kxe t& ) d x F^>{k,x 



J8\ 



Hi 







(14) 



which makes d x $£ continuous even if the potential has a sharp cut-off at a point x = xq. 

Substituting this ansatz, Eq. flT2|), into Eq. (^) and using Eqs. ( |T3"D and (fH|), lead to 
coupled differential equations 



d x Ft\k,x, 



e -.(-) 



— ti[\kxe ld ) V(xe ie ) 
2ik 



hf\kxe i0 ) F^ , (k,x,e) + h [ f'(kxe t0 ) Ff\k,x 



(-)/ 



(15) 



and 



d x Ft\k,x, 



e -,(+), 



-—hY\kxe ie ) V(xe ie ) 
2ik 



x 



hf\kxe w ) F^ > (k,x,9) + h [ f , (kxe td ) Ff\k t x 



J6\ 



A6\ 77l(— ) , 



with boundary conditions 



Ft\k,0,6) = Ff ) (k,0 



(-)/ 



(16) 



(17) 



that follow immediately from Eqs. (|H]) and (TTj. With these boundary conditions, this 
system of coupled differential equations are equivalent to the pair of integral Volterra-type 
equations, 



F^\k,x,9) = 1 + 



J 8 



ik Jo 



h^ikx'e^Vix'e^eik, x'e i0 ) dx' 



(18) 



and 



Ft\k,x,9) = 1 - ~— I hf } (kx'e lti )V{x'e w )<5> e {k,x'e w ) dx' . 



ik Jo 



(19) 



Iterated solution of these equations, on a small interval [0,x m j n ] and with ji(kxe ie ) as 
the zero-order approximation to &i(k,xe t6 ), can be used to obtain the boundary values 
F^ ± '(k,x min ,9), which in turn are needed to solve the differential equations, Eqs. ( |i~5"D and 
(|j~6|), numerically. Thus, we can find the complex dilation $^(/c,xe ) of the regular solution 
in the form Eq. (|12D , using either differential (Eqs. flTop and (0)) or integral (Eqs. ([18]) 
and (|TP|)) equations. 

The question then arises as to the physical meaning of this solution and that of the 
functions F^ ± \k,x,9). To answer, we note that the complex rotation, Eq. (0), makes the 
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Hamiltonian nonhermitian, so it can have complex eigenvalues ||10|| . Indeed, the Jost function 
can have zeros, not only on the imaginary axis of the (complex) fc-plane but also in the entire 
lower half of that plane. Schematically this is depicted in Fig. 1. Therein, the points on the 
positive imaginary axis (fc i, k 02 , k 03 ) represent Jost function zeros corresponding to physical 
bound states of the system. The points shown in the fourth quadrant (^04,^05); are the 
resonance state zeros. The dividing line defined by the rotating angle 9 is also displayed. 
This line is the limit for the resonance position angle p for which the condition Im kr > 
is fulfilled. 

The asymptotic form of the function &e(k , xe l9 ), corresponding to these zeros, decreases 
exponentially if the angle 9 is chosen large enough, i.e. greater than the angle (p of Eq. (^j), 
since || 

hf\kxe ie ) ^=Fiexp[±i(fcce* - y)]. (20) 

Hence <&e(ko, xe %e ) are square-integrable eigenfunctions of the rotated Hamiltonian corre- 
sponding to either real (bound) (Ihko = 0,lmk > 0) or complex (Ftk > 0,lmk < 0) 
eigenvalues. Although these eigenvalues, being independent of 9, always correspond to the 
energies of the physical bound and resonance states, the corresponding solutions §i(k , xe l0 ) 
acquire physical meaning only when 9 = 0. However, if one is interested only in eigenen- 
ergies, it is convenient to use 9 > thereby dealing with square-integrable wave functions. 
This idea was successfully used for the location of resonances with variational methods |J. 

In the present approach, the square-integrable nature of the dilated regular solution at 
resonance energies is not used directly. Instead, the crucial feature for our approach is that 
the function h^^xe 10 ) decreases exponentially at large x if Im(kxe l6 ) > 0. This is true 
for the whole region of the complex fc-plane in which the bound and resonance state zeros 
lie (see Fig. 1). In contrast, the regular solution is square integrable only at the discrete 
points k 0i of this region. According to Eq. ([TT]) , at all other points k not coinciding with the 
spectrum k 0i of the rotated Hamiltonian, the regular solutions diverge as 

<S> e (k,xe i0 )~ e \ Im ( kxei(h >\ . (21) 



As can be seen from Eqs. fll8Q and (19), both functions at the spectral points fcoi, have 
finite limits when x — > 00 



\imF( + \k 0i7 x,9) = ft(k* oi ) , (22) 



and 



limF, ( - ) (A; 0l ,x,^) = /,(^)=0, (23) 



which must coincide with the asymptotic coefficients defined in Eq. ([TT|). At any other point 
k ytz k Qi within the zone Im (kxe 19 ) > 0, the integral of Eq. ([19|) still has a finite limit, 



\im Ft\k,x,9) = f e (k) , (24) 

because hf + \kxe l6 ) is a decreasing function. However, the integral in Eq. (|18|) diverges. 
This is not surprising and by no means a defect to our approach. Simply, it reflects the 
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well-known fact (see, for example, Refs. and [11]) that the Jost solutions, T\ (k,r), of 
the Schrodinger equation with boundary conditions, 



lim e Tikr rt ±] (k,r) = 1 , (25) 

r— >oo 

are defined within different regions of the complex fc-plane : Im k < for T% and Jm A; > 
for t\ . The only common region of their definition is the real axis Imk = i.e. the 
physical region. Only for a special class of rapidly decreasing potentials it would be a band 
along this axis 0,p|ll|. However here we consider the most general case defined by the 
condition Eq. (fj). 

The same is true for the limit (x — > oo) of our functions Ff ~°\k, x, 8) with 6 = 0, which 
are closely related to the Jost solutions, as one can see from Eq. flT2|). However, in contrast 
to the Jost solutions, the functions F^(k, x, 0) are well defined for any finite x, because the 
boundary conditions, Eq. (|l7l), are imposed on them at the point x = 0. 

When 9 > 0, the border dividing the two regions of complex k where the limits of 
F^ ± \k, x, 6) exist, is rotated downwards to the line \k\e~ ld . The fact that above this line the 
function F^ + \k, x, 9) has no limit for x — > oo does not contradict the asymptotic requirement 
of Eq. flTT|). Written in general form, these asymptotics must include both terms, 

h[ + \kxe ie )ft(k*) and h [ f ] } {kxe ie )f £ {k), 

only on the dividing line, Im (kxe l9 ) = 0, where they are of the same order of magnitude. 
Above this line, Eq. ( |TTD is a sum of 'small' and 'large' terms; the first of which can be 
omitted. The 'small' term, h^\kxe ie )f^(k*), becomes significant above the dividing line 
only at the discrete points k = k 0i , i = 1, 2, . . . , N where the 'large' term, ^^(kxe 19 ) fe(k) , 
disappears because fi(koi) = 0. And, as we have already seen, at these points the limit 
Eq. ( p2]) exists and the 'small' term is well defined. At all other points the growth of 
Ff + \k,x,8) is suppressed by the exponentially decreasing factor h^\kxe ld ) and the 'small' 
term becomes negligible as compared to the 'large' term. Indeed, at a sufficiently large x 
the integrand of Eq. (|Tj]) behaves like 

~ exp{-2ikxe i9 )V{xe i9 ). 

If the potential at infinity is a nonzero constant, then the integration would produce F^ ~ 
exp(— 2ikxe 10 ). However, the potential obeys the condition Eq. (§) so that F^ diverges 
less rapidly than exp(— 2ipxe ld ). Therefore the product h^F^ grows less rapidly than 
^^(kxe 10 ) and does not affect the asymptotic properties specified by Eq. ([□]). The two- 
term decomposition, Eq. (|T2|), guarantees term by term correspondence with the asymptotic 
form Eq. (|ll]) at the discrete points k$i and also on the dividing line Im (kxe 10 ) = 0. For the 
other points (k ^ koi, Im (kxe 10 ) > 0), the first term of Eq. (|l2j) contains some admixture 
df (k,x,6) of the diverging term. However, this admixture is infinitesimal as compared to 
the leading term, h[ \kxe' l9 )fe(k) of the long range asymptotic behaviour, 

<f>t{k, x, 6) ^ ~ {[h { +\kxe ie )fZ(k*) + dt(k, x, &)} + h { f\kxe te )f,(k)} , (26) 
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i.e. 

Jim {[h^ikxe^f^k*) +d e (k,x,6)} / (kxe l9 )f e {k)} = . (27) 

The asymptotic behaviour of the functions F^ + \k, x, 6) and F^~\k, x, 6) below the dividng 
line | A: | e — 461 is just the opposite, i.e. is finite and ^ diverges. However, the function 
Fg \k,x,8) has a physical content only when Im(kxe l6 ) > 0, while Fg (koi,x,9) at the 
discrete points k 0i} i = 1, 2, . . . , N, where F^~\k i, oo, 9) = 0. 

Thus, if the potential permits a large enough choice for 6, we can calculate the Jost 
function for all complex values of k ^ of interest, including the momenta on the positive 
imaginary axis and in the resonance region under the positive real axis as well. As a conse- 
quence of Eq. (|24D , the Jost function can be found by solving the coupled one-dimensional 
equations, Eqs. (0) and (|16|), from x = to x = x max , where x max is determined by the 
required accuracy and on how fast the potential decreases. The faster it decreases the sooner 
d x Fg \k,x,6) vanishes and F^ \k,x,8) tends smoothly to the constant fe(k). 

It is interesting to note that at sufficiently large x the system of equations, Eqs. (|15D 
and (|16D , decouple as the 'small' term h^F^ becomes insignificant as compared to the 
'large' term h^Fjf . 

III. NONCENTRAL SHORT-RANGE POTENTIAL 

In this section we consider a two-body problem with spin-dependent interaction which 
does not conserve the orbital angular momentum I. As in the previous section, we assume 
that the matrix elements 

V(jMe's'),(jMes)( r ) =< yis' I yfs > 

of the potential operator, sandwitched between the spin-angular momentum eigenfunctions 

mil 

of the total angular momentum operator, 

j = e + s, 

obey the same constraints (0) and (|2|) at the origin and infinity. 

To simplify the notation, the superscripts JM will be omitted and instead of the pair of 
subscripts is, satisfying the triangle condition 

\J - s \ < £ < \J + s\, 

we will use the symbol [£}. 
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A. Boundary conditions 



The physical wave-function, describing the state with a given momentum k, is the product 

r) = -(y. u ) = -J2 yi M (r)u [e] (k, r) 
r 

of the [£]-row of y and the [£]-column u of the radial components satisfying the coupled 
differential equations 

D e (k,r)u[£\(k,r) = ^ V[i> W] {r)u[ V ]{k, r) (28) 

m 

with (physical) boundary conditions 

0) = 0, 

(29) 

where the asymptotic function U[e](k,r) has a different form for the bound, scattering, and 
resonance states. 

As it is more convenient to deal with universal boundary conditions imposed at a single 
point, we consider the so-called fundamental matrix, r), of the regular solutions 

instead of the physical column solution, U[^(k,r). The physical solution then can be con- 
structed as a linear combination of columns of (A; , r). By definition, $ is a fundamental 
matrix if it has the following three properties ]2[]: 

i) Each column of is a solution of Eq. (|28|), i.e 

D t (k,r) $[mt>](k,r) = J2 V mn( r ) ®{t"W]( k i r ) > ( 30 ) 

[£"] 

ii) (A;, 0) = 0, and 



iii) all columns of &[e][£] are linearly independent. 
Moreover, must be a square matrix since Eq. ( p8|) has as many independent regular 

solutions as the column dimension [12 . 

The property ii) implies that the desired universal boundary condition should be imposed 
at the point r = 0. We emphasize here that the boundary condition $[i'][g\ (k, 0) = 0, is not 
sufficient to define this matrix since Eq. (|30D is of second order. Furthermore, Eq. ([Tj|) is 



singular at the origin, and consequently the existence and uniqueness theorem pi is not 
valid at this point. However, this theorem is valid for any small r > 5 > 0. Hence, the 
matrix $ can be uniquely defined by 8) and d r <&(k,5), while within the interval [0,6], 
$(k,r) can be evaluated explicitly as follows. After multiplying by r 2 and using Eq. (|J), 
Eqs. (P0|) decouple giving 



r%-£(£+l) $ [m (k,r) = 0, re [0,5] 



(31) 
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Each of these equations has two independent solutions, namely the functions ~ r e+1 and 
~ r . Therefore, near the origin, we can use the Riccati-Bessel function, 

*w ^ 2^|) r ' +1 • < 32 > 

as the regular solution. Hence, as a first order approximation to $\g\\m(k,r) for r e [0,5] we 
can take the diagonal matrix 

®lW']( k > r ) r ~ Q S [W'} 3t'(kr) , (33) 

satisfying the above-mentioned three properties. Higher order corrections give non-zero 
off-diagonal elements, so to guarantee linear independence of the columns, the off-diagonal 
elements of a row must be infinitesimal in comparison to the diagonal element of that row. 
The boundary condition thus reads 

$ m >](k,r) _ 
r->o ]£'{kr) 

which is a generalization of Eq. (|5]). 

Since $ is the fundamental matrix, any physical solution is a linear combination of its 
columns with r-independent coefficients, i.e. 

u [e] (k,r) =^<$> m , ] (k J r)A [e] (k) . (35) 

[t>] 

The coefficients Am must be chosen to satisfy the appropriate boundary condition Q2§| ) at 
infinity. 



B. Complex rotation 

To generalize the approach described in the previous section, we perform the complex 
rotation, Eq. fl7|), of the coordinate. The corresponding complex dilation of the fundamental 
matrix, $[^/](A;, xe ld ), is the solution of 



D t {ke i9 ,x) $ [m (k,xe w ) = e™Y,V m ^xe w ) Qp^k, xe») 

[£"] 



(36) 



where 



lim 



$ [m (kxe i0 ) lj t »(kxe™) 



i0\ 



(37) 



We then apply the variable-constant method i.e. look for a solution of the form 

$ [m (k,xe i9 ) = X - [h { +\kxe ld ) + h[-\kxe w ) F^(k,x,e) 



(38) 



with 
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h { +\kxe ie ) d x F^, } (k, x, 9) + h['\kxe ie ) d x F^ £l] (k, x,0) = 0. 



(39) 



The resulting coupled equations for the new unknown functions, Fr^? e Jk, x, 9), have the form 



,W 



^h[-\kxe*)Y,V mn 



xe 



i0\ 



i6\ 



and 



id 

' hi + \kxe* e )J2V mn (xe w ) 



(40) 



2ik 

(+) 



hg/t^(kxc ^)-^\£"i\i r \ (^"j ^) ~\~ hpn ^ (kxc )-F 1 j£//^^ (fc, x, ^) 



(41) 

The appropriate boundary conditions for them can be derived from Eqs. (|37|) and fl38|) . First 
though, one must compensate the singularities of h^ikxe 10 ) and h[~\kxe l6 ) at x = 0. For 
this we can demand F^ and F^ to be identical as x — > 0, viz. 



Ffifa (k, x, 0) F [em (k, x, 9) 



and 



as then 



Ffafa (k, x, 9) (A;, x, 9) 



$ m >](k, xe w ) x ~ Q 3e(kxe ty ) F [m (k, x, 9) 
so that, from Eq. fl37j), the boundary conditions become 



lim 



j£(fcxe^)F^] ;] (fc, x, #) / >(fcxe 



S0\ 



(42) 



(43) 



(44) 



(45) 



However, as x = is a singular point, for practical calculations one needs to solve Eqs. (]40[) 
and ( |4l| ) analytically on a small interval (0, 5} and then impose the boundary conditions at 
x = 5. Such an analytical solution can be found easily by choosing 5 to be small enough so 
that for x G (0, 5} we can write 



d x F<&(k,x,e) = ±—h [ t'(kxe M ) V [m {xe w ) Mkxe") . 

IK 

For small x, in h t = j£ ± m^, the Riccati- Neumann function is dominant, so that 

^^(fc, x, 9) ~ -— n,(fcre*) %*](ze*) >(fcxe ie ) . (47) 
Upon integration of this (approximate) equation one finds, 



-,i0 



(46) 
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F^k, x, 9) = - e — J n e {kxe ie ) V [m {xe ie ) j £ >(kxe ie ) dx + const , (48) 

and if the arbitrary constant of integration is taken to be equal to this function has 

the asymptotic behaviour, 

Ffg^k, x, 9) x ~ o 5 m >] --jTJ Mkxe ie ) V [m {xe ie ) j v {kxe i6 ) dx , (49) 

which obeys the condition of Eq. (f£5|). In practical calculations the last indefinite integral 
can be found analytically by using the leading terms of series expansions of ri£, Vj^/j, and jp. 
In this way one can find that the second term of Eq. (^) is regular at x = if t < £' (right 
upper corner and the diagonal of the matrix) and may be singular if i > £' (left lower corner 
of the matrix). But this singularity is compensated by the presence of je in Eq. (f44|). 

Thus, the coupled equations, Eqs. ( |40"D and fl4T|), along with the boundary conditions, 
Eqs. (|45| ) and (f49f ), form a well-defined differential problem. 

By analogy with the case of a central potential, we find that F^? £ Jk, x, 9) has a finite 
^-independent limit, 

\imF^ f] (k,x,e) = f^ f] (k) , (50) 

when Im (kxe t6 ) < (i.e. in the complex fc-plane below the dividing line), while F^n has 
the limit, 



lim F, 



We 



^(k,x,9) — /[£][]/] (k) , (51) 



when Im(kxe ie ) > (i.e. in the complex /c-plane above that line). Therefore, the funda- 
mental matrix has the asymptotic form, 

$ m] (k, xe* e ) — > \ [hf\kxe i6 ) f^ f] (k) + h[-\kxe w ) , (52) 

with fw^n{k) being the matrix generalisation of the Jost function. We refer to this as the 
Jost matrix. 

Substituting Eq. (|52|) into Eq. (|35|), one finds that the physical solutions involve only the 
first term of Eq. (52) if 

J2f^ f] (k)A m (k) = 0. (53) 

[£'] 

This system of homogeneous algebraic equations defining the coefficients A\g\ has a nontrivial 
solution if and only if 

det||/^(/c )|| = 0. (54) 

Therefore the zeros k Qi of the determinant of the Jost matrix are the spectral points corre- 
sponding to the bound and resonance states. The square-integrable complex-dilated wave 
function of any such state is 

*(k ,r,9) = %L £ yjMQ) [h^\k Q xe ie ) F£&(ftb,M) 



2r 

AX [£][£'] 



+h ( f\k xe w ) F^(ko,x,6)\ A m (k ) . (55) 
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IV. POTENTIAL WITH COULOMBIC TAIL 



Consider the case where the leading term of the long-range behaviour of the potential 
has the form 2-qk/r. We split the potential as 



V{r) = {V(r) - ^} + 2*2 



(56) 



so that the entry in brackets is a short-range interaction. We designate that hereafter as 
V sr {r) and note that it satisfies the constraints Eqs. ([I]) and (§). If the radial Schrodinger 
equation now is written as 



D](k,r)<&t(k,r) = V sr (r)<f> e (k,r) 



(57) 



where 



D](k, r) = d 2 r + k 2 - 1(1 + l)/r 2 - 2r]k/r , (58) 
the natural generalisation of the boundary condition, Eq. (H), for the regular solution is 

\im[$ e (k,r)/F i (r ] ,kr)} = l , (59) 

with Fi(r), kr) being the standard regular solution of the Coulomb equation §|, 

D](k,r)F e (r],kr) = . (60) 
Application of the complex rotation, Eq. (0), gives 

D](ke t6 , x) $ e (k, xe ie ) = e 2ie V sr (xe ie ) <S> e (k, xe w ) , (61) 

and 



lim 



$i(k,xe ld ) /F e (r), kxe 



1 



(62) 



Then, following the procedure developed in Sec. II, solutions $>e(k, xe zd ), are sought that 
have the form 



&e(k, xe 



i0\ 



1 r 



H { e + \ri,kxe te ) F^ + \k,x,9) + H i i -\r ] ,kxe ie ) F^~\k,x,9)} , 

where now the Lagrange condition is 

H { e +) {7], kxe id ) d x Ff + \k, x, 6) + H^(r], kxe w ) d x Ff~\k, x, B) = . 
Therein are the combinations of the regular and irregular Coulomb functions, 

HP{r),z)=F e (r],z)TiG £ {r),z) , 
having the asymptotic behaviour, 



(63) 



(64) 



(65) 
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H^\tj,z) — ► =p exp(±i{z — n In 2z — — + argT(£ + 1 + irj)\ 



\z\— »oo 

In the limit 77 — ^ 0, these are reduced to the Riccati-Hankel functions h[ (z). 

The above process gives first order coupled equations for the F^ <y± \k, x, 9), namely 



(66) 



d x F? i+ \k,x,9) = ^-Ht\v,kxe* 9 ) V sr (xe w ) 



H ( e + \ri,kxe w )F; i ^ ) (k,x,6) + Hf > (7 1 ,kxe w )F; i( ~>(k,x,6) , (67) 



J0\ j?v(+) 



rH/ 



i6\ r?v{-), 



and 



d x F^~\k,x,9) 



Hi +) {r],kxe ty ) V sr {xe w ) 



2ik 



x 



H { e + \r],kxe M ) F^>(k,x,9) + Hf>(r ] ,kxe iy ) F r e l{ -'(k,x, 



J9\ ttV(+) 



(-)/ 



,i0\ TpV(-), 



(68) 



It is easily seen that the asymptotic properties of F\ (k,x,9) are the same as those of 
F^ ± \k,x,9) discussed in Sec. II. Therefore the limit 



lim F^\k,x,9) = p t {k) , 



(69) 



exists for all complex k obeying the condition Im (kxe ) > (i.e. above the dividing line). 
The Coulomb modified Jost function is thus fi(k), and it follows that the corresponding 
long range limit of Fg'(k, x, 9) exists in the zone Im (kxe l6 ) < and for all spectral points 
k 0i where f^(k 0i ) = 0, % = 1, 2, . . . . 



V. A NUMERICAL EXAMPLE 

To illustrate the method, we consider s-wave solutions for the exponential potential, 

2VTl 

V(r) = T vexp(-r/R) . (70) 

For this case, the differential problem, defined by Eqs. (|]) and (H), can be solved analytically 
0. The corresponding Jost function is given by 

2ikR 



/„(*)= T(l - 2ikR) J. 2lkR ( 2R ] j \ , (71) 

and is defined for all complex k except at the discrete points k n = —in/ (2R), n = 1, 2, 3, . . . 
on the imaginary axis. 

There are two parameter sets for this potential that are of particular interest. The first, 
with strength v = 155.17MeV and range R = 0.76 fm, describes the triplet iViV-interaction 
13] ( /i 2 /(2/i)=41.47 MeV-fm 2 ). With this interaction, the Jost function has a zero on 



the positive imaginary axis in the complex A;-plane, corresponding to the deuteron binding 
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energy. The second parametrisation, is v = 104.2 MeV and R = 0.73 fm for the singlet NN- 
interaction \ 13^ . In this case there is a zero of the Jost function on the negative imaginary 
axis, corresponding to virtual singlet deuteron. More specifically, these zeros of the function 
(0) lie at k = i 0.23640 fm" 1 (-2.3175 MeV) and at k = -i 0.039913 fm' 1 (-0.066064 MeV) 
respectively. The relevant solutions, Fq (A;, x max , 9), of Eqs. (|15|) and (fL6|), obtained by 
means of numerical Runge-Kutta integration up to x max = 20 fm, have zeros at the same 
momenta for all choices of 9 e [0, ~). 

The successful evaluation of the virtual state via Fq (A;, x max , 9), despite the fact that 
the corresponding zero lies below the dividing line, is due to the short range nature of the 
potential. In fact, for such potentials the non-dilated Jost solutions are both well defined 
within a narrow band along the dividing line 0; a band wide enough to include the zero 
linked to the virtual state. 

To check the equivalence of Fo~\k, x max , 9) to the exact Jost function for different values 
of complex momenta and to observe ^-independence of Fq \k,x max ,9) above the dividing 
line, we used the triplet deuteron potential (first set of parameter values) and obtained the 
solutions for six values of the momenta that are situated on the arc of the radius \k\ = 2 fm -1 . 
With k = \k\ exp(iip n ) , we chose (p n to be it/2, 7r/4, 0, — vr/8, — vr/4, and — 37r/8. The 
calculations were repeated with four different values of the dilation angle 9 (0, tt/8, 7r/4, 
and 37r/8 ) which correspond to different directions of the dividing line (see Fig. 2). A point 
\k\ exp(i(f n ) is above this line when Im (kxe ) > 0, i.e. when ip n + 9 > 0. 

The results are presented in Table I. For non-negative values of 0, the calculated functions 
agree very well (to 1 part in 10 5 ) with the exact Jost function and the solutions are very 
evidently ^-independent. For negative <ft however, close agreement with the exact Jost 
function values and ^-independence of the calculations only occurs with the choice 9 > \<p\. 
As expected from our theory, the divergence is extreme when this condition is violated. 

These results demonstrate the effect of the complex rotation in that the greater one 
can take 9, the wider is the sector of the fourth quadrant of complex /c-plane where 
Fq (A;, x max , 9) coincides with the (exact in the test case) Jost function and therein the 
calculated function is independent of 9. 

In Table II, the results of our calculations of the test example are shown for different 
coordinate l x' values in the range 0.5 < x < 20 fm but at a fixed momentum and with 9 
set to zero. These results clearly display the divergence of F (+) (A;,a;,0) and the convergence 
of Fq \k,x,0) when x — > oo. The third column of this table shows a rapid convergence 
of Fq ^ which reflects the fact that this function is smoother than the whole function $o 
and varies only at distances where the potential is significant. Indeed, Fq~^ has converged 
very well ( 1 part in 10 5 ) by x = 10. But the relevant asymptotic terms are the ratios 
of the terms in Eq. (|i"2|) . As i — > oo, these ratios decrease rapidly. The last column 
of Table II explicitly shows how the first ('small') term, h^" Fq, of the asymptotics of 
Eq. (P7|), becomes infinitesimal in comparison with the second ('large') term, /i^Fo - -*, 
despite the divergence of Fq + ' (k,x,0). 
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VI. CONCLUSIONS 



We have shown how a combination of the variable-constant and complex coordinate 
rotation methods can be used to recast the two-body Schrodinger equation with the diverse 
boundary conditions to consider bound, virtual, scattering, and resonance states, into a set 
of linear first order coupled equations for Jost-type solutions, with which all possible state 
forms can be treated in a unified way. The derived equations are amenable to numerical 
intergration to obtain the Jost functions (Jost matrices) for all momenta of physical interest. 
The viability of this approach was demonstrated by its application to find solutions for a 
potential for which the (s-wave) Jost function has an analytic form. 

An advantage of the proposed systems of equations, over studies of two-body states in 
which the Schrodinger or Lippmann-Schwinger equations are used, is that the bound state 
and scattering problems can be treated by the same method. In addition, resonance states 
can also be handled. Moreover, for the scattering problem a much higher accuracy can be 
achieved because both and are smooth functions since the oscillatory behaviour 
of the asymptotic form of the complete solution $ is factored out in this approach. 

With our approach, one can examine the contribution of different parts of the potential 
curve in the scattering process, simply by truncating the integrations at different distances. 
That is equivalent to truncation of the potential. This information is contained in F^' 
( as a function of r ) because it is equal to the Jost function for the truncated potential at 
distance r. 

Another important advantage of our approach over the conventional ones, is that one 
can calculate the Jost function for complex momenta in the resonance region. To our 
knowledge, this is the simplest way to locate resonances (zeros of det F^~'(k, x max , 9)), with 
high accuracy. Such information is of crucial importance in various fields of Physics. Inverse 
scattering problems where the parametrization of the S matrix is required, is a case. In the 
^-dependent Marchenko inversion method ||, for example, the S matrix is written as 

S , W = ^ = n |±^, (72) 
fi(k) V k-a n 

and the interpretation of the parameters and their location in the complex fc-plane are 
important in order to obtain the correct potential. 

The proposed method has certain drawbacks however. One is that when 9 ^ 0, analytic 
continuation of the potential to the first quadrant of the complex r-plane is required. That 
may not be easy to do when the potential function of interest has discontinuities or is given 
in the form of a numerical table. A second drawback concerns the point k = at which 
the proposed equations are singular. This problem can be overcome by using a prescription 
given in Ref. |TJ]]. Within a small region around k — 0, the Riccati-Hankel functions h\ (At) 
can be expanded in a power series || and each term therein factorized in k and r. Similarly 
can be expanded in powers of k with unknown r-dependent coefficients in this region; 
the coefficients being specified by the resulting system of coupled, k- independent, differential 
equations. 

The range of application of the proposed method can be extended further. Since the 
derived formulae are valid for non-integer values of the angular momentum £, the method can 
be used to solve the iV-body hyperradial equation of the hyperspherical harmonics approach 
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15] and which has the form of the radial Schrodinger equation but with half integer values 



for £. Using non-integer £, we can also consider potentials that tend as ~ r~ 2 toward the 
origin, by including the singularity of the potential into the centrifugal term. Such cases 
arise when supersymmetric (SUSY) transformations are applied to remove a bound state 
from the spectrum. The resulting supersymmetric partner potential has a ~ 1/r 2 behaviour 



at the origin |T6| 



Finally, we can consider complex values of the angular momentum £, and, in principle, 
study all the problems concerning Regge trajectories. 
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TABLES 



TABLE I. Comparison of Fq (k,x = 20fm.,#) with the exact Jost function (Eq. (|7l|)) at the 
six momentum points situated on the arc defined by k = 2 e itp fm _1 . The numbers given in the 
brackets give the exponent power of ten. 





Hp) 












<p 


e = o 


9 = - 

8 




e = — 

8 


exact 


TT 


Re 


0.56706 


0.56706 


0.56706 


0.56706 


0.56706 


o 
z 


Im 


U. 


— U.oz { — o) 


n i a i <\\ 

u.(4 [— yj 


U.lz 4J 


U. 


TT 


Re 


0.59205 


0.59205 


0.59205 


0.59205 


0.59205 


1 

4 


Im 


— u.ziybo 


— U-ziyoo 


—u.ziybo 


n ii n£Q 
— U.ziybo 


—u.ziybo 


n 


Re 


0.69491 


0.69491 


0.69491 


0.69490 


0.69491 


u 


Im 


-0.48975 


-0.48975 


-0.48975 


-0.48972 


-0.48975 


TT 


Re 


12.3 


0.81233 


0.81233 


0.81231 


0.81233 


~8 


Im 


-7.5 


-0.66767 


-0.66767 


-0.66764 


-0.66767 


TT 


Re 


-2 (11) 


-87 


1.0360 


1.0360 


1.0360 


~4 


Im 


3 (12) 


47 


-0.85549 


-0.85546 


-0.85549 


3vr 


Re 


9 (17) 


3 (12) 


1 (4) 


1.5184 


1.5185 


~~8~ 


Im 


9 (19) 


-2 (13) 


3 (4) 


-1.1802 


-1.1802 



TABLE II. Values of Fq and F for different x. The values of k and 6 are fixed at 
2 e J 4 fm. _1 and respectively. The numbers in brackets give the exponent power of ten. 



(fm) 


Ft\k,x,Q) 


r o 


(k,x,Q) 




0.5 


0.62 + i 0.23 


0.82292 


- i 0.06269 


0.20 


1 


0.59 + z 1.40 


0.70057 


-i 0.15146 


0.13 


2 


3.30 - i 2.24 


0.62144 


-i 0.20521 


0.02 


3 


-12.4 + i 11.6 


0.59993 


- i 0.21596 


0.6 (-2) 


5 


-353 + i 23.8 


0.59261 


- i 0.21939 


0.4 (-3) 


10 


5 (4) + i 7 (5) 


0.59205 


- i 0.21965 


0.6 (-6) 


15 


(9) - i (8) 


0.59205 


- j 0.21965 


0.8 (-9) 


20 


-2 (11) - i 3 (12) 


0.59205 


- i 0.21965 


0.1 (-11) 
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FIGURES 

Im k 




FIG. 1. Schematical picture of a typical distribution of the spectral points koi on the complex 
A;-plane. The dividing line defined by the rotation angle 9 is also shown. 
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Im k 




FIG. 2. Several dividing lines corresponding to the Jost function calculations given in Table 
I. The points at which this function was calculated are shown by the dots. 
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